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Abstract. We calculate the beaming of coherent emission 
by means of geometrical optics. When spherical elemen¬ 
tary waves are phase-coupled in that way that constructive 
interference occurs along one axis, the intensity is lower for 
off-axis viewing angles. The spatial angle of constructive 
interference shrinks to a fraction of 47 t when the length 
scale of the emission region is larger than the wavelength 
of the radiation. This leads to an amplification of the re¬ 
ceived radiation in addition to the relativistic beaming 
caused by the relativistic radial outflow which is thought 
to be important for the radio emission of pulsars. The ef¬ 
fect will be calculated by numerical methods and can be 
approximated by a simple analytic expression. 
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The main feature of a coherent process is that there 
are N particles radiating in phase (e.g. Melrose 1991 and 
references therein). Therefore, the coherent intensity of 
the radiation is calculated as N 2 times the single particle 
amplitude. In the case of N incoherently radiating parti¬ 
cles the interference terms cancel and the total intensity is 
given by N times the square of the single particle ampli¬ 
tude. Consequently the intensity of N coherently radiating 
particles is N times higher than in the incoherent case. 

But, as will be shown, if the volume of coherently radi¬ 
ating particles (hereafter coherence volume) is about the 
same size or larger than the wavelength of the emitted 
radiation the phase coupling needed for the increased in¬ 
tensities can only be maintained in one direction, i.e. along 
one axis. Seen from off-axis angles there may even be 
destructive interference leading to zero intensity. In this 
contribution we calculate the angle at which the intensity 
drops to zero. Thus we get a reasonable estimate for the 
beaming due to coherent radiation. 
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Assuming that the reason for coherence is a free elec¬ 
tron maser (FEM) process (i.e. inverse Compton scat¬ 
tering (ICS) of relativistic electrons and nonlinear elec¬ 
trostatic plasma waves (Asseo et al. 1990; Benford and 
Weatherall 1992; Asseo 1993)) we derive a formula for the 
beaming angle depending only on the density, the Lorentz 
factor of the electrons and the size of the coherence vol¬ 
ume. 

Consider a volume of N particles which are able to 
radiate in phase. For simplicity we assume that the vol¬ 
ume is a cube with one of its axes in the direction where 
fully constructive interference occurs. This seems reason¬ 
able since nonlinear electrostatic plasma waves appear 
pancake-shaped due to Lorentz contraction (along the 
magnetic field lines) and the electrons are moving per¬ 
pendicularly to the extension of these plasma solitons. 

Therefore it seems natural to expect that the direction 
of relative velocity between waves and particles describes 
the direction of phase coupling (this shall be chosen to be 
the z- axis being parallel to the vector (0,0,1)). The prob¬ 
lem can be reduced to a two-dimensional one by orientat¬ 
ing the coordinate system in such a way that the observer 
is located in the direction (sin a, 0, cos a) and the y- axis 
can be disregarded in the calculations. 

If the coherence volume is d 3 (where d is the length 
of a cube edge) we can choose the origin of the coordi¬ 
nate system to be a corner such that in the cube all three 
coordinates are normalized to the range of [0..d]. 

Throughout this paper we perform our calculations in 
the beam frame, i.e. a Lorentz factor belongs to the non¬ 
linear wave structure approaching the electron beam (ap¬ 
proximately equal to the Lorentz factor of the beam seen 
by an observer). 

The term coherence volume needs some clarification. 
It is the volume occupied by the coherently radiating par¬ 
ticles which is the same as the volume of the considered 
soliton. In other words, we do not distinguish between the 
part of the soliton involved in the reaction and the vol¬ 
ume of beam electrons radiating coherently except for the 
Lorentz contraction since the volume is observed in the 
beam frame. The soliton can be treated as the cause of 
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coherence whereas the beam particles are coherently radi¬ 
ating when they interact with the soliton. 

Next we calculate the phase shift of the wave from a 
radiating particle having the coordinates (x, z ) relative to 
a similar particle at the origin (0,0) when radiating at an 
angle of a. (The y-coordinate is disregarded due to the 
choice of the coordinates). The phase difference can be 
expressed as 

<j>(x, z) = k (S(x, z) + z) (1) 

where k = 2 tt/\ is the wave number of the emitted ra¬ 
diation, 8{x,z) denotes the difference in pathlengths for 
radiation from the two particles (where the usual approx¬ 
imations are adopted for calculations of interference at 
optical grids). The second term inside the brackets, z, de¬ 
scribes the intrinsic phase difference which is necessary for 
constructive interference in z-direction. Obviously a par¬ 
ticle which is right ahead of another in the distance of z 
must have a phase delay of kz to interfere constructively. 



Fig. 2. The same as Fig. 1, but for two particles having the 
same x- coordinate. 



Fig. 1. Schematic picture showing how the path difference for 
two points is calculated having the same z- coordinate. 


d 0 is now calculated for two special cases: 

1) z=0: It can be seen from Fig. 1 that do = a: sin a ss 
xa for small values of a. 

2) x=0: In this case we find do = —2 cos a (cf. Fig. 2) 
and thus </>(0, z) = z (1 — cos a) « za 2 /2 for small a 

We add these two phase differences to get the general 
result. This renders: 


<j>{x,z) 


27T 

T 


xa + - za“ 


( 2 ) 


the sum by an integral does not affect the results heavily 
when TV 1. Using the integral has the advantage that 
particles need not to be arranged in a grid (as the elec¬ 
tron temperature might be quite high) and the result will 
become more representative of the physics involved. 

It is reasonable to assume that coherence is not pro¬ 
vided instaneously at its maximum possible number N, 
but that the coherence volume will be zero initially and 
then expand in some way. Since we are mainly interested 
in the maximum possible effect, we do not assume a certain 
time dependence but just set a maximum size of the co¬ 
herence cell. The coherence is caused by a plasma process, 
therefore the most natural assumption is to set the size to 
d = c/u pe . Here v pe = \J(ne 2 ) / (m e £oj) / (2n) is the (frame 
invariant) plasma frequency (cf. Weatherall and Benford, 
1991). 

By re-scaling the integration interval to [0..1] (that 
means transforming the integration variables to dimen¬ 
sionless lengths in units of the coherence cell extension) 
we finally obtain for the relative phase averaged wave am¬ 
plitudes in the direction of a 


A 

A) 



dz da; 


( 3 ) 


where A is the amplitude of the wave and A 0 = A(a = 0) 
the corresponding amplitude of waves emitted parallel to 
the particle propagation. 

The double integral simplifies to 


For calculating the total intensity relative to the max¬ 
imum value we integrate the amplitudes over one phase 
and over the coherence volume and take the square of 
the integral. In principle one has to sum up the complex 
amplitudes of waves from all the particles but replacing 


A 

Ao 


2 

M 2 a 3 


(sin a sin 6 


cos a cos b + cos a + cos b — 1)(4) 


where M = 2ir c/ (Xv pe ), a = Ma, b = Ma 2 /2. The first 
zero of the above equation yields the angle of the minimum 
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and we find the solution a = w — b which leads to a simple 
quadratic equation for a. The final result then reads 


C^min 





(5) 


The last approximation is valid as long as M 1 which is 
always the case as we consider only plasma processes for 
radio emission, and thus can write c/(Aiy pe ) = v em /v pe = 

where g is an emission model dependent parameter 
which is at least 1. 

This result can be compared to a more exact geometric 
model that we obtain by numerical integration over the 
phases of waves emitted from a soliton. 

The soliton shape is chosen since a coherent plasma 
process requires efficient bunching of the radiating parti¬ 
cles. Physically speaking this is the same as a high density 
modulation or a strong longitudinal plasma wave, which 
is well described by an almost non-dispersive longitudinal 
plasma soliton with a certain extension perpendicular to 
the magnetic field lines. 

We assume a cylindrically symmetric shape like 
cosh (-()) “ (see Fig. 3) where p denotes the distance from 
the center of the soliton in units of its characteristic length 
(called ’’extension” of the soliton in the simplified discus¬ 
sion before) and R as the characteristic length scale equal 
to c/v pe or equal to d, used in the simple approximation 
above (e.g. Weatherall 1997). 

The solitons developing from electrostatic plasma 
waves are negligibly thin in the direction of the field lines 
but have a large perpendicular extension. Let the observer 
be placed at a large distance s and at an angle a with re¬ 
spect to the normal direction. Again we assume that the 
positive x- axis is the projection of the observer’s line of 
sight. Then, s being the distance and cf> the meridional 
angle in the soliton plane we find the expression 


r 2 = s 2 — 2sp sin (f> cos a + p 2 . 


( 6 ) 


A spherical wave emitted from a point source can be writ¬ 
ten as 

p -ikr 

vk(r) = — (7) 

Thus the superposition of all elementary waves from a 
soliton reads 


oo 2tt 

C r p -i kr(s,p,ot t <f>) / O \ ~ 2 

A so i(s, 9) = A 0 / / —-—pcoshf —J dip dp. (8) 

J J r [s, P, ot, <p) ^ 77 ' 

0 0 

k is the wave number with the same definition as 
above. This equation can be treated numerically. We find 
that the result is not altered significantly but only slightly 
modified by the different shape of the coherence cell. As 



Fig. 3. Typical shape of a soliton. The maximum amplitude 
of the soliton is normalized to 1, and points in radial direction 
(z). The plane spanned by the line of sight and the direction 
of the outflowing particles as assumed to be the y=0 plane 
. X- and y- axis are scaled in units of the characteristic length 
(see text). 


can be seen in Fig. (4) the beaming is even stronger than 
found in the simple approximation. 1 

As mentioned above we assume a plasma process to 
be the origin of coherent emission. As an illustrative ex¬ 
ample we take the free electron maser mechanism which 
can (cum grano salis) be treated as inverse Compton scat¬ 
tering of plasma solitons and relativistic electrons. Then 
the incoming electron sees a plasma wave with a wave¬ 
length that is equal to that measured in the soliton frame 
(which is almost the same as the rest frame for very low- 

1 No great discrepancy arises even though we treated the 

’’coherence cell” as a cubic box in the first case and in the 
numerical calculations we assumed the soliton to be pancake¬ 
shaped. But if the phase-coupled emitter extends over many 
wavelegths in all directions, \/d will dominate the interference 
effects and the symmetry of the emitter will be of secondary 
importance. 
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Intensity versus angle for gamma=5 



Fig. 4. Beaming pattern for the coherent emission. The exact 
result is shown by the fully drawn line whereas the dotted 
line shows the approximation. The angle between the line of 
sight and the direction of movement is given in degrees and 
the intensity is normalized to 1 for a = 0 


relativistic solitons) divided by 7 due to relativistic space 
contraction. Seen in the beam frame the electrostatic wave 
is simply reflected with the same frequency. In the rest 
frame the reflected wave has a frequency that is higher 
by another factor of 7 and is given by 7 2 a> pe as can be 
verified by a Lorentz transformation from the beam frame 
to the observer’s frame. A wave of this frequency cannot 
propagate as a plasma wave any more and will be decou¬ 
pled from the plasma propagating as a free electromag¬ 
netic wave. Because we performed our calculations in the 
beam frame, only one 7 is required in the formula for the 
emitted frequency; that means <7 = 1 in the formula given 
above. 

Thus we finally find the expression 

Ctmin 71 (9) 

27 

From that equation we obtain that the beam becomes 
sharper when 7 is larger. Because we require the strongest 
coherence at low frequencies (and therefore naturally low 
Lorentz factors) it is natural to assume that we have 
non-stationary processes which start incoherently (at high 
Lorentz factors) and become more and more coherent 
when 7 decreases. Equation (9) should only be valid for 


the lowest Lorentz factors when coherence is fully devel¬ 
oped. 

The angle of the minimum now provides an estimate 
for the spatial beaming angle due to coherence. This is 
given by 

AH = na 2 min (10) 


which is true for small angles. This corresponds to a rela¬ 
tive beaming factor of 


4w is 2 

m = An = lfh 


(ii) 


in the beam frame. 

For a fixed observer in addition to coherent beaming 
there is, of course, beaming by a factor of 47 2 / 7 t (Ry- 
bicki and Lightman 1979) due to the relativistic motion 
of the solitons themselves. Thus an observer sees a relative 
beaming factor that reads 


47T 9 

^ An 7 


64 7 4 
7r 


( 12 ) 


It has been shown above that in coherent processes 
one can find a natural beaming effect that is caused by 
the non-isotropic emission pattern of a phase-coupled grid 
of point sources due to interference. This leads to higher 
fluxes of coherent radiation in special directions which can 
be important for estimations of brightness temperatures 
or observed radiation ” power” (that means the power that 
would have to be emitted if the radiation was isotropic 
in the observers frame, see Manchester and Taylor 1977). 
The described effect might be a contributional cause of 
the very intense fluxes of radio photons in pulsars. Some 
detailled discussion of the flux problem will be subject of 
a later paper. 
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